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ABSTRACT 



Context. Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously 
reported for M31 X-ray sources, observed by XMM-Newton. However, we have recently discovered that these PDS 
result from the improper addition/subtraction of non-simultaneous lightcurves. 
Aims. To understand the properties and origins of the artefact. 

Methods. We re-analysed our XMM-Newton observations of M31 with non-simultaneous and simultaneous lightcurves, 
then combined simulated lightcurves at various intensities with various offsets and found that the artefact is more 
dependent on the offset than the intensity. 

Results. The lightcurves produced by the XMM-Newton Science Analysis Software (SAS) are non-synchronised by 
default. This affects not only the combination of lightcurves from the three EPIC detectors (MOS1, MOS2 and pn), 
but also background subtraction in the same CCD. It is therefore imperative that all SAS-generated lightcurves are 
synchronised by time filtering, even if the whole observation is to be used. We also find that the reported timing 
behaviour for NGC 4559 ULX-7 was also contaminated by the artefact; there is no significant variability in the correctly- 
combined lightcurves of NGC 4559 ULX-7. Hence, the classification of this source as an intermediate-mass black hole 
is no longer justified. 

Conclusions. While previous timing results from M31 have been proven wrong, and also the broken power law PDS in 
NGC 4559 ULX-7, XMM-Newton was able to detect aperiodic variability in just 3 ks of observations of NGC 5408 
ULX1. Hence XMM-Newton remains a viable tool for analysing variability in extra-galactic X-ray sources. 
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1. Introduction 

The variability and spectral properties of Galactic low mass 
X-ray binaries (LMXBs) are well known to depend more on 
the accretion rate than on whether the pr imary is a neu- 
tron star or black hole (Ivan derKlis 1994). At low accre- 
tion rates, the power density spectra (PDS) may be char- 
acterised by broken power laws, with the spectral index, 7, 
changing from ~ to ~1 at some break frequency in the 
range 0.01-1 Hz (jvan der Klislll994fk such variabi lity has a 
r.m.s power of ~10-40% ( e.g. Ivan der K lis 1995). We de- 
scribe such PDS as Type A ((Barnard et al.ll2004f ) . At higher 
accretion rates, the PDS is described by a simple power law 
with 7 ~1 and r.m.s. variability <10% (van der Klis 1994 
1995h . We refer to these PDS as Type B ((Barnard et all 



2004). 



Type A variability is characteristic of disc-accreting 
X-ray binaries; any X-ray source that exhibits such vari- 
ability may be identified as an X-ray binary, rather than 
a foreground star or background active galaxy (see e.g. 
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iBarnard et al.|[2003bD . Such variability has been reported 
in XMM-Newt on observations of M31 by s ome of the au- 
thors (see e.g. IBarnard et all l2003bl 12003: IWilliams et ail 
120051 ). Type A PDS have also been report ed for ultralu- 
mino us X-ray sources (ULXs) in NGC 4 559 (|Cropper et al.l 
[2001 and NGC 5408 dSoria et al.ll200l . 



However, we have now discovered an artefact that pro- 
duces false Type A PDS, caused by the improper addi- 
tion of lightcurves from the three EPIC instruments on 
board XMM-Newton (MOS1, MOS2 and pn), where the 
lightcurves are not synchronised (Barnard et al., this is- 
sue). The data were analysed using SAS 6.0.0 and FTOOLS 
5.3.1. In the Sect, d] we discuss how the error occurred. We 
then describe a new, corrected analysis of the M31 data in 
Sect.O Section [4] describes our investigation into the arte- 
fact using simulated lightcurves. Then, Sect. 5 details the 
causes of the artefact. Section [5] re-analyses XMM-Newton 
observations of two ultraluminous X-ray sources with pub- 
lished Type A PDS. We draw our conclusions in Sect. [7] 
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Fig. 1. Power density spectra of 0.3-10 keV lightcurves of r3-60, using Method I (left) and Method II (right). We present 
best fits to the Method I PDS of power law and broken power law models; a broken power law is clearly required. The 
PDS is hence reminiscent of Type A variability in Galactic LMXBs. However, the Method II PDS is flat, revealing that 
this variability is artificial. 



2. Manipulating lightcurves generated by SAS 

XMM-Newton has o b served M3 1 extensively, and 
iBarnard et all (|2003fJbl I2004J 120061) conducted a sur- 
vey of bright sources in the central region, over four 
observations. For each observation of every source, 0.3-10 
keV lightcurves were obtained from source and background 
regions in the MOS1, MOS2 and pn images, using 2.6 
second binning. 

We have recently discovered that the results of 
manipulating lightcurves produced with the SAS tool 
evselect using the lcmath FTOOL produces differ- 
ent results, depending on how the time selection in 
evselect is phrased. Originally, lightcurves were pro- 
duced from events files that were filtered in evselect 
with the expression "(PI in [300:10000]) && (TIME 
in [tstartdstopf' ; additional filtering was instrument de- 
pendent: MOS lightcurves were additionally filtered us- 
ing a #XMMEA_EM&&(PATTERN<=12}\ while the pn 
equivalent was "#XMMEA_EP&&(PATTERN<=4)" . We 
shall refer to this filtering as Method I. 

However, different lightcurves are obtained when evs- 
elect is used in either of the following ways. Firstly, one 
may assign tstart and tstop to the TLMIN1 and TLMAX1 
keywords in the pn, MOS1 and MOS2 events files, us- 
ing FTOOLs such as fmodhead or fparkey; we refer to 
this method as Method II. Method III involves filtering 
each lightcurve with the evselect tool, using the parame- 
ters "timemin= tstart timemax= tstop" . Methods II and III 
are equivalent. Vitally, they are not equivalent to Method 
I, although this is not mentioned in any documentation. 
Consequently, we analysed lightcurves obtained with each 
method to determine which, if any, were correct. 

3. Re-analysing the data 

3.1. Combined EPIC lightcurves 

We obtained MOS1, MOS2 and pn lightcurves from the 
source extraction regions of several sources in the M31 ob- 
servations, using Methods I and II. For Method I, we found 
that uncertainties of combined EPIC (MOSl+MOS2+pn) 



lightcurves were underestimated by a factor of ~10-40%. 
This ratio varied within each lightcurve, with the standard 
deviation in the ratio for a given lightcurve decreasing with 
increasing luminosity. The reason for the small uncertain- 
ties is discussed in Sect. 5.1. 

Closer inspection of the lightcurves has revealed the dif- 
ference between Method I and Method II lightcurves: the 
Method II lightcurves are synchronised, while the Method 
I lightcurves are not. Of course, for the Method II or III 
lightcurves to be synchronised, tstart must be the same for 
all detectors, and must be present in all the events lists. 

Many PDS from lightcurves combined using Method I 
exhibit Type A variability, while identical data combined 
with Method II do not. To investigate the cause of the Type 
A variability in the observed PDS, we examined the Method 
I and II lightcur ves of XMMU J004208.9+412048 (r3-60 
iKong et al.il2002D . The Method I and Method II PDS are 
shown in the left and right panels respectively of Fig.[IJ they 
were obtained using the FTOOL powspec. The Method I 
PDS is well modelled (% 2 /dof = 6/8) by a broken power 
law where the spectral index changes from to 0.33±0.05 
at a frequency of 29±4 mHz. Such a PDS is characteristic 
of Galactic LMXBs. However, the PDS of the Method II 
lightcurve is flat, with a Leahy power of 2, consistent with 
Poisson noise. Hence the variability in the Method I PDS 
is artificial. 



3.2. Background Subtraction 

Background subtraction requires two lightcurves: source 
and background. We produced background-subtracted pn 
lightcurves of various X-ray sources from the 2002 January 
XMM-Newton observation of the M31 core, to test whether 
the lightcurves from the same detector were synchronised. 
We present the Method I PDS from source, background 
and background subtracted pn lightcurves of r3-125 in 
Fig- 121 panels (a), (b) and (c) respectively. The source and 
background PDS are flat, while the background-subtracted 
PDS is of Type A. Hence, background subtraction causes 
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Fig. 2. PDS from source+background (a), background (b) 
and background-subtracted (c) pn Method I lightcurves of 
r3-125. The background-subtracted PDS is of Type A, al- 
though the component PDS are flat. 

artefacts when non-synchronised source and background 
lightcurves are used. 

4. Investigating the artificial variability with 
simulated lightcurves 

To determine the cause of artificial variability we com- 
bined artificial lightcurves with varying intensities and 
offsets. We defined 100 evenly spaced intensity intervals 
from 0.01 to 1.00 count s -1 , and for each interval gener- 
ated a single "source" lightcurve and 100 "background" 
lightcurves (using the terminology of lcmath: the "back- 
ground" lightcurve may be added to or subtracted from the 
"source" lightcurve). This intensity range reproduces well 
the observed intensities from the bright X-ray sources in the 
M31 Core. Each lightcurve had 60 ks lightcurve duration, 
with 2.6 s binning. 

For each lightcurve, we generated a series of photons 
with random arrival times, using the Park & Miller Minimal 
Standard random number generator, shuffled by the Bays & 
Durum method (Knuth, 1981); the random number seed for 
the nth "background" lightcurve at intensity i was defined 
as — (lOOOOi + n). The source lightcurves were generated 
with seeds of — lOOOOi; hence, each lightcurve was repro- 
ducible and unique. We then conducted a series of tests for 
artificial power in the PDS of combined lightcurves. 

4.1. Adding two lightcurves with different start times 

Firstly we added each "background" lightcurve for every 
intensity interval to the corresponding "source" lightcurve 
at various offsets with lcmath; the offset was defined as the 
delay in seconds as a fraction of the 2.6 s time binning, and 
ranged from 0.05 to 0.95 at intervals of 0.05. We then pro- 
duced geometrically grouped, 64-bin PDS of the resulting 
lightcurve with the FTOOL powspec, averaged over ~180 
time intervals. Initially, we examined the PDS for excess 
variability by fitting them with a constant Leahy power of 
2; for our PDS, a x 2 >25 / 10 dof signified a 3er detection 
of power above the Poisson noise. We find that the offset 



is more important than intensity, with all combinations of 
lightcurves with offsets 0.25-0.80 producing significant ex- 
cess variability; 12% of lightcurves with an offset of 0.05 
or 0.95 produced artificially variable PDS. However, not all 
variable PDS are considered Type A. 

To be classed as Type A, a power law fit to the PDS 
must be rejected, and a broken power law fit must be rea- 
sonable. A broken power law fit may be described by spec- 
tral index a that changes to spectral index (3 at frequencies 
higher than some break frequency the final parameter is 
the power at iv When fitting the PDS, we rejected power 
law fits with good fit probability <0.005 (\ 2 >24) and ac- 
cepted broken power law fits with good fit probability > 
0.05 (x 2 < 16); for the broken power law fit, we set a=0, as 
observed in Galactic LMXBs. Less than 10% of lightcurves 
with offsets in the ranges 0.05-0.25 and 0.75-0.95 exhibited 
Type A variability, while 30% of lightcurves with / = 0.50 
were of Type A. 

4.2. Adding three lightcurves with different start times 

Since Barnard et al. combined source lightcurves from the 
three EPIC detectors in the original M31 survey work (e.g. 
iBarnard et aIll2003rJlal[200l we added a third lightcurve to 
the combined artificial lightcurves discussed in the previ- 
ous section. For Case A, this lightcurve was the 50th back- 
ground lightcurve with an offset of 0.1, and for Case B 
the 50th background lightcurve with an offset of 0.5 was 
added. The resulting lightcurves had intensities of 0.03-3 
count s -1 , with offsets of 0.05-0.95 as before. 

We tested the resulting PDS for Type A variability. The 
results are summarised in Table[TJ Again, we found that the 
fraction of Type A lightcurves dependent on offset rather 
than intensity, with the maximum at / = 0.50. However, 
even an offset of 0.05 produced artificial Type A variability 
in 12% of combinations for Case B. Combining three non- 
synchronous lightcurves (in Case B) results in a factor of 
>2 more false Type A variability than just combining two 
lightcurves. Hence the artefact is strengthened as more non- 
synchronous lightcurves are added (combining detectors) or 
subtracted (background subtraction). 



5. The causes of the artefact 

Three factors contribute to the artificial Type A variability: 
underestimated uncertainties, excess power in the PDS, and 
suppression of high-frequency variability. We now discuss 
these in turn. 



Table 1. Summary of variability observed when three arti- 
ficial lightcurves are combined, for Cases A and B. For each 
case, the fraction of combined lightcurves exhibiting Type 
A were shown, for offsets / = 0.05, / = 0.50 and / = 0.95. 



Case 


/ = 0.05 


/ = 0.50 


/ = 0.95 


A 


0.003 


0.13 


0.011 


B 


0.12 


0.75 


0.12 
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5.1. Underestimating Method I uncertainties 

The smaller uncertainties of Method I are due to the work- 
ings of lcmath. Whether adding or subtracting lightcurves, 
lcmath treats the first lightcurve as the source and the 
second lightcurve as the background; it calculates the 
mean "background" intensity for each time interval in the 
"source" lightcurve. If the lightcurves are offset by some 
non-integral number of time bins, then the intensities of 
two "background" intervals are weighted by their overlap 
with the "source" interval. The intensity of the n th interval 
of the summed lightcurve, S (n), is then given by 



S (n) = A (n) + / x B (n - 1) + (1 - /) x B (n) 



(1) 



where A is the "source" lightcurve, B is the "background" 
lightcurve and / is the fractional offset. The latter is given 

by 



' n 



AT 



(2) 



where t$ and tg are the start times for the source and 
background lightcurves, respectively, and AT is the width 
of the time interval. 

By default, the lightcurves begin at the arrival time for 
their first photons. In practise, however, t$ is determined by 
the arrival of the first photon in the "source" lightcurve but 
lcmath sets tg to the start of the first "background" inter- 
val to overlap the "source" lightcurve. In addition, lcmath 
ignores all source intervals with no "background" interval. 
Hence, < / < 1. For synchronised lightcurves, t$ and tg 
are both set to the same user-defined value. Hence / = 
and S (n) = A (n) + B (n), as it should be. 

The underestimated uncertainties are caused by an er- 
ror in the lcmath code lcmathrdbkg.f. When two or 
more background bins coincide with a source interval, with 
count rate C\ and fractional overlap Oi, then lcmathrd- 
bkg.f assigns an uncertainty of 0\ x Cf , rather than the 
Gaussian uncertainty (Oi x C;) ' 5 if no uncertainty is pro- 
vided. Similarly, when an uncertainty o\ is provided, then 
lcmath calculates the total uncertainty as [(Oi xo- 1 ) 2 +(0 2 x 
cr 2 ) 2 + - ■ -] a5 , when it should be (O x x erf + 2 x <j\ + . . .) ' 5 . 
Since Oi < 1, Oi/O® 5 <1, and the uncertainties are under- 
estimated. 

The underestimated Method I uncertainties are also re- 
sponsible for the excess power in the PDS. This is because 
the Leahy PDS produced by powspec are normalised such 
that the white noise level expected from the data errors cor- 
responds to a power of 2 (see powspec documentation). To 
do this, Gaussian errors are assumed and the normalisation, 



.4 



Leahy 3 



must be 



- 9 A * 



(3) 



where At is the sampling time, N is the number of points 
and cr 2 is the mean squared error (see e.g. Appendix A of 
IVaughan et al.|[2003t ) . I.e., A Lea hy oc 1/er 2 . Therefore, under- 
estimating the Method I uncertainties by a factor of 1.1—1.4 
will overestimate the PDS power by a factor of ^1.2-2. 

5.2. Suppression of high-frequency power 

After correcting for the power excess described above, we 
find that the Method I PDS actually dips below the Poisson 



noise level at high frequencies. I.e., the high frequency vari- 
ability is suppressed, causing the break in the PDS at fre- 
quency V\,. 

For those M31 sources identified with broken power 
law PDS, v h = 20±4-57.8±0.5 mHz, independent of the 
luminosity. We have studied the break frequencies of 
the 190,000 lightcurves that were constructed from three 
non-simultaneous, randomly-generated lightcurves and dis- 
cussed in Sect. 14.21 We found to vary over ~3-70 mHz, 
with no strong correlation between and offset or inten- 
sity. While the break frequency is necessarily related to the 
bin size, this is relation is very weak also because Vh varies 
by a factor of ~30 while the bin size is constant. 



6. Ultra-luminous X-ray sources in NGC 4559 and 
NGC 5408 

We conducted a brief literature search for published PDS 
of combined EPIC lightcurves that may also be affected by 
this data- analysis artefact. We found a published broken 
power law PDS for NGC 4559 U LX7 (ICropper et al.l l200l 
and also for a ULX in NGC 5408 (|Soria et al.ll2004h . We re- 
analysed XMM-Newton observations of these sources, first 
following the methods described in these papers, using non- 
synchronised lightcurves, then comparing their PDS with 
those from synchronised lightcurves. However, it is beyond 
the scope of this paper to re-analyse all published XMM- 
Ne wton variability data ! 

ICropper et alj (|2004l ) reported a broken power law PDS 
from an ultra-luminous X-ray source (ULX) in NGC 4559 



with a UV/X-ray luminosity exceeding 2.1x10 erg s 
We extracted MOS1, MOS2 and pn lightcurves from a 
source region with 30" radius using Method II. We con- 
structed a PDS that was normalised to give fractional 
r.m.s. 2 variability, and the expected noise was subtracted. 
The expected noise had a power of 4.3. The frequency bins 
were geometrically grouped. The resulting PDS is shown 
in Fig. [3] This PDS is consistent with zero power (good 
fit probability = 0.22); hence, the variability reported by 
Cropper et al. ( 2004) is artificia l. 

We note that lFeng fe Kaaretl (|2005| ) also re- analysed the 
XMM-Newton observation of NGC 4559 U LX- 7, and were 
able t o reproduce the variability reported by I Cropper et ail 
(|2004l ) only when the intervals of background flaring were 
included; hence, they mistakenly attribute d the variabil- 
i ty to the background flares, unaware that ICropper et al.l 
(|2004[) also excluded t he bac kground flaring from their anal- 
vsis. iFeng fe Kaaretl (|2005f ) produced their PDS directly 
from the events lists that describe each photon (Kaaret 
2006, private communication). Hence they did not observe 
the artefact after the intervals of high background were ex- 
cluded 

ISoria et all (|2004h studied a series of XMM-Newton ob- 
servations of a ULX in the dwarf galaxy NGC 5408 over 
an 18 month period between 2001, July and 2003, January. 
Its 0.2-12 keV spectrum requires two components, a power 
law with photon index T = 2.6-2.9 and a thermal compo- 
nent with blackbody temperature of 0.12-0.14 keV; such 
a spectrum is typical of a black hole X-ray binary in the 
very high / steep power law state (|McClintock fc Remillardl 
2004 ^. The inferre d 0.2-12 keV luminosity is ~10 40 erg s" 1 



( Soria et alj |2004). They present a broken power law PDS 
from ~3 ks of the 2003, January XMM-Newton observation 
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of NGC 5408, where 7 increases from -0 to 1.3±0.2 at 2.5 
mHz. Such a PDS is also consistent with a black hole X-ray 
binary in the very high state. 

We extracted EPIC li ghtcurves from a circular source 
region with 30", lollowing lSoria et all (|2004h . A Method II 
PDS is presented in Fig. [4j the PDS is r.m.s 2 normalised, 
and averaged over 3 intervals of 256 bins, with 5 s time bins. 
Fitting this PDS with zero power results in a x 2 /dof of 
92/14, hence the variability is real. Figure [4] also shows the 
best power law and broken power law fits to the PDS, with 
X 2 /dof = 77/12 and 24/11 respectively. The broken power 
law fit is clearly preferred (f-testing yields a 99.95% chance 
that the improvement is significant), but we caution that 
it is not an acceptable fit: the goodness of fit probability is 
just 1.2%. 



brighter sources will naturally have smaller offsets between 
CCDs. 

The published broken power law PDS of NGC 4559 
ULX7 is also artificial. The PDS of the NGC 5408 
ULX shows genuine intrinsic variability. Timing analy- 
sis of XMM-Newton observations of extra-galactic X-ray 
sources is indee d viable. M31 X-ray s ources have exhib- 
ited pulsations (lOsborne et al.l | 2001t iTrudolvubov et all 
120051) : bursts (jPietsch k. Haberll 12003): periodic dipping 
due to photoelectric absorption (jTrudolyubov et all 12002: 



Mangano et al 



([Barnard et al 



200 



2006 



from a processing disc in one case 
finally branch movement in a tri- 



modal c olour-intensity diagra m reminiscent of Galactic Z- 
sources (jBarnard et al.ll2003aD . 



7. Conclusions 

Some of the authors have previously reported variability in 
XMM-Newton observations of M31 X-ray sources that are 
characterised by broken power law P DS, which are signa- 
tures of disc-accreting X-ra y binaries (jBarnard et al.l l2003b 
12004 IWilliams et all 12005)'. However, this "Type A" vari- 



ability observed in XMM-Newton observations of M31 is an 
artefact and may be entirely attributed to errors introduced 
by lcmath when combining non-synchronised lightcurves. 
(Barnard et al., this issue). However, extragalactic ana- 
logues of the Galactic X-ray binaries should exhibit Type 
A and Type B behaviour; Type A variability should be ob- 
served in sufficiently bright X-ray sources. 

In the survey of XMM-Newton observations of X-ray 
sources in M31, we found that the Type A variability was 
observed at lower luminosities than Type B variability, as 
expected for LMXBs. However, the most likely explanation 
for this observed behaviour is that the offset was smaller 
for the brighter sources; the start time for an XMM-Newton 
lightcurve is defined by the arrival time of the first event, so 
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